Photoionization of Fe XV 
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Relativistic and resonance effects in the photoionization of Mg-like Fe XV are investigated using 
the Breit-Pauli R-matrix Method (BPRM) at near-threshold and intermediate energies, comple- 
mented by the Relativistic Random Phase Approximation (RRPA) and multi-channel quantum 
defect theory in the energy region up to the L-shell ionization thresholds. The cross sections exhibit 

. . . extensive resonance structures that considerably enhance the effective photoionization of Fe XV. 

0^ , These results should be of general interest in photoionization modeling of X-ray sources observed 

O^ ■ by space observatories. 
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There is particular interest in the photoionization of 
highly charged iron ions owing to their importance in the 
f~^ ■ modeling of astrophysical plasmas such as in active galac- 
P\J I tic nuclei and quasars, supernovae, stellar coronae etc., 
and in laboratory sources such as inertial and magnetic 
confinement fusion devices [1]. Among the most com- 
plex atomic species are those isoelectronic with the third 
row elements containing open n = 3 shells. The K- and 
Xf-\ ] L-shell photoexcitations and photoionizations are highly 
^^ ■ topical areas of investigation as they are generally in the 
X-ray range accessible to satellite observations, such as 
On from the recently launched Chandra X-ray Observatory 
^^ (CXO), and the upcoming X-ray Multi- Mirror Mission 

^ , (XMM) and Astro-E [1]. 
Qh' While previous calculations exist for the iron ions, they 

J,^ ' are in relatively simpler approximations that neglect the 
5^ ■ multiplet or the fine structure, and autoionizing reso- 
"ti ' nances that are known to enhance the cross sections sig- 
C^ nificantly [2]. Recently, extensive photoionization calcu- 
►^ lations were carried out in the close coupling approxi- 
• i-H . mation using the R-matrix method under the Opacity 
^^ ■ Project [0P,3] for most astrophysically abundant ele- 
H , ments. While the OP cross sections include the resonance 
■ - - ' effects, they are in LS coupling neglecting fine structure 
effects in the detailed cross sections, as shown in this 
work through relativistic close coupling calculations us- 
ing the Breit-Pauli R-matrix codes developed for the Iron 
Project [4]. 

Whereas the low energy region can be accurately con- 
sidered with the BPRM method, at higher energies lead- 
ing up to the L- and K-shell ionization the number of 
coupled channels becomes very large and the computa- 
tions quite prohibitive. On the other hand, there are 
very few studies done of the resonance structures lead- 
ing up to these ionization edges [5]. It is therefore of 
interest to not only carry out the low energy BPRM cal- 
culations, but also to study the nature of resonances at 
high energies using other methods such as the relativis- 
tic random phase approximation (RRPA) and relativis- 
tic multi-channel quantum defect theory (RMQDT) . The 



aim is to combine the two techniques to investigate pho- 
toionization cross sections over the entire energy range 
of practical interest, as a guide to more elaborate close 
coupling calculations. 

The computationally intensive BPRM calculations are 
carried out as described in recent works [6] . The residual 
(target) ion following photoionization is represented by 
an TV-electron system, and the total wavefunction of the 
initial bound (A^+l) electron- ion system is represented 
in terms of the ion cigcnfunctions as: 

^{E)=AY,X^O^+Y.''J^l^ (0-1) 

i 3 

where Xi is the ion wavefunction in a specific state 
SiLiJiTTi and 9i is the wavefunction for the (7V-|-l)-th 
electron in a channel labeled as SiLi{Ji)7rik'^£i [Jtt]; 
kf being its incident kinetic energy. The ^j's are the 
correlation functions of the (Af+l)-clcctron system that 
account for short range correlation and the orthogonal- 
ity between the continuum and the bound orbitals. The 
<I>j's may also give rise to bound channel resonances at 
intermediate energies due to inner-shell photoexcitation- 
autoionization, as described in this work. 

The eigenfunctions Xi oi the residual ion Fe XVI are 
obtained through atomic structure calculations using the 
program SUPERSTRUCTURE [7]. Table 1 hsts the Xi 
corresponding to 17 levels dominated by configurations 
2p^ {^S) n£, up to n = 5. The calculated eigenenergies 
are compared with the cxpcrimetal values. Although the 
Fe XVI levels included in the close coupling expansion 
span the n = 5 levels, most of the relevant resonance 
structure appears below the first excited fine structure 
levels ^Pol/2, 3/2. However, owing to the high ion charge 
the n-complexes are well separated and the coupling ef- 
fects, including resonances at higher energies, are weak. 

We consider the photoionization of the ground level of 
Fe XV, 3s2 (i^o) + hiy — >e + Ss (^Si/2), coupled to the 
continua of the ground and the excited levels of Fe XVI. 
Fig. 1 shows the detailed cross section with the Ryd- 
berg series of resonances in the region up to, and slightly 
above, the first two excited thresholds 3p (^Pi/2 3/2)- ^'^^ 
comparison, the non-rclativistic cross sections from the 
Opacity Project work by Butler et al. [8] arc also shown 
(dashed lines). It is seen that the OP results, in LS cou- 
pling, did not resolve the extensive resonance structures 
obtained in the present calculations including relativis- 
tic fine structure. In the region ^P°i2 ~^ ^zn ^^^ ''^S" 
onances are considerably weaker than in the region be- 
low the ^Pii2 owing to autoionization into the excited 
3p ^P°i2 continuum, in addition to the ground level con- 
tinuum 3s ^5'i/2)- Fig. 2 shows the photoionization cross 
section in an extended energy range that spans all ex- 
cited Fe XVI levels in the close coupling expansion (Table 
1). Above the 3p ^P°i2 3/2 thresholds however, resonance 
structures are fairly isolated and relatively far apart. The 
background cross sections and some resonance positions 
agree well with the earlier LS coupling OP results [8]. 
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FIG. 1. Photoionization cross section of the ground level of 
Fe XV 3s (5*0) in the near-threshold region; solid line - rel- 
ativistic BPRM results, dashed line - non-relativistic results 
[OP,8]. 
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FIG. 2. Comparison of the present BPRM cross sections 
with the earlier OP non-relativistic calculations over an ex- 
tended range spanning all target thresholds in the close cou- 
pling expansion. 

BPRM calculations can be further extended in the in- 
termediate energy range, from the highest target level 
to the n = 2 inner-shell ionization thresholds. This re- 
gion does not include the target levels that would cor- 
respond to the first sum in the close coupling expansion 
(Eq. [1]). However, a number of (N-l-l)-electron config- 
urations, ^j{Fe XV), are included in the BPRM calcu- 
lations in the second sum in Eq. [1]) with eigenvalues 
in this range. These Fe XV configurations correspond 
to outer 3s electron excitation, as well as the inner 2p 
electron excitations, i.e. [ls^2s^2p^] 3i^ + 3i, nt, n < 
5,£ < 2, and the inner-shell excitation configurations: 
[ls^2s^2p^] U'^nt,n < 5,£ < 2. Although the 2p-shell 
ionization are not explicitly represented in the BPRM 
calculations, the 2p-shell excitation- autoionization is in- 
cluded to some extent via resonances corresponding to 
these ^j{Fe XV). Some of the resonances in the 30 - 80 
Ry region, due to inner-shell photoexcitations, are much 
stronger than those due to the valence 3s-shcll excitations 
in Figs. 1 and 2, as discussed below (Fig. 4). 

Close to the ionization threshold of the n = 2, we 
carry out RRPA-I-RMQDT calculations, with the de- 
tailed resonance strucutures shown in Fig. 3. The 
RRPA and RMQDT calculations are carried out as 
described in [9] and [10]. In the autoionization res- 
onance spectrum region between 80-86 Ry, the pho- 
toionization cross section in open dissociation channels, 
that is, the 3s photoionization cross section, is calcu- 
lated and plotted as a function of photon energy in 
Fig. 3. The seven interacting relativistic dipole chan- 
nels considered for the RRPA-I-RMQDT calculations are: 
3s — > e(p3/2),e(pi/2);2p3/2 — > nd^ / 2, 3 / 2, ns 1/2; and 
2pi/2 — ^ '^c?3/2j^si/2- The theoretical single configu- 
ration Dirac-Fock thresholds for 3s,2p^/2,2pi/2 are at 
33.314, 86.26, and 87.21 Ry respectively. Below the 2^3/2 
threshold there are five interacting Rydberg series cor- 
responding to discrete excitation from 2^3/2 and 2pi/2 
levels. Coupling between these excitation channels and 
the ionization channels from the 3s levels causes config- 
uration interaction between the discrete and continuum 
states, which in turn leads to autoionization resonances 
seen in Fig. 3. 
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FIG. 3. RRPA and RMQDT cross sections with resonances 
converging on to the L-shell ionization thresholds. 



Finally, Fig. 4 shows the complete set of BPRM and 
RMQDT+RRPA results obtained. The strongest res- 
onances in the middle range between 54 - 70 Ry cor- 
respond mainly to inner-shell excitation configurations 
2p^ 3s^ {3p,3d), particularly to the 2p — > 3d exci- 
tation. We identify these resonances from an exten- 
sive atomic structure configuration-interaction calcula- 
tion for Fe XV using SUPERSTRUCTURE, employ- 
ing the same one-electron orbital basis set as the tar- 
get ion Fe XVI. For example, the autoionizing lev- 
els 2p5 3s2 3p 3'i(S'i,Po,i,2, -00,1,2,3), lie between 54.49 
- 56.83 Ry. The large resonance complex between 
58.26 - 59.76 Ry corresponds to the autoionizing lev- 
els 2p5 3s2 M ^•H'S'i,Po,i,2,-Do,i,2,3)- The experimen- 
tal first ionization potential of Fe XV is 33.58 Ry (the 
present BPRM calculated value is 33.56 Ry), and the 
highest Fe XVI target threshold is at 25.41 Ry. This im- 
plies that the L-shell resonance complex just below 60 
Ry coincides with the highest target levels 5^(^1)3/2,5/2) 
in the wavefunction expansion (Table 1). It is clear that 
the inner-shell photoexcitation-autoionization would be a 
major contributor to the effective photoionization cross 
sections below the L-shell ionization threshold(s). 



TABLE I. Fe XVI target level energies (Ry) in the 
wavefunction expansion 



Index 


Configuration [Term] 


J 


E(Obs) 


E(Calc) 


1 


2p''{'S)3s ['S] 


1/2 


0.00000 


0.00000 


2 


2/(iS)3p pP°] 


1/2 


2.52596 


2.52186 


3 


2/(i5)3p pP°] 


3/2 


2.71688 


2.69095 


4 


2p'^CS)3d fD] 


3/2 


6.15544 


6.21637 


5 


2/(i5')3rf pD] 


5/2 


6.18198 


6.25459 


6 


2/(iS)4s p5] 


1/2 


17.0182 


17.0836 


7 


2/(i5)4p pP°] 


1/2 


18.0252 


18.0643 


8 


2/(i5)4p pP°] 


3/2 


18.0980 


18.1269 


9 


2/(i5')4d pL»] 


3/2 


19.3570 


19.4013 


10 


2p^CS)4d [^D] 


5/2 


19.3677 


19.4186 


11 


2p«(i5')4/ pF°] 


5/2 


19.9077 


19.9638 


12 


2pS(i5')4/ pF°] 


7/2 


19.9125 


19.9703 


13 


2/(15)55 p^] 


1/2 


24.2500 


24.3164 


14 


2/(i5)5p pP°] 


1/2 


24.7606 


24.8020 


15 


2p''CS)5p pP°] 


3/2 


24.7970 


24.8357 


16 


2p«('S)5d pL>] 


3/2 


25.4065 


25.4552 


17 


2p«(i5')4d pD] 


5/2 


25.4116 


25.4639 
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FIG. 4. The BPRM and (RRPA + RMQDT) cross sec- 
tions in the energy range up to L-shell ionization. The arrow 
at 59.98 Ry is the highest target threshold in the BPRM cal- 
culations. The total sub-shell photoionization cross sections 
a(2p) and a{2s) are computed using the RRPA. 



The RRPA and RMQDT results in Figs. 3 and 4 show 
an excellent match with the background BPRM cross 
sections. The resonances just below the 2p ionization 
threshold form a pseudocontinuum whose effective area 
corresponds to the total photoabsorption in this region, 
which in turn rises to the threshold value at the 2p ioniza- 
tion energy of 86.26 Ry. The total photoabsorption cross 
section should be continuous across ionization thresholds. 
Therefore the cross section at and above thrcshold(s) is, 
in principle, comparable in magnitude to the effective 
resonance averaged photoexcitaion-autoionization cross 
section below threshold(s). The energy region up to the 
2p and 2s ionization thresholds in Fig. 4 is dominated by 
resonance structures, many of which are not yet consid- 
ered. The inner-shell ionization edges in particular are 
not sharp jumps, as often obtained in non-resonant cal- 
culations, but as diffuse pseudocontinua rising up to the 
relevant threshold cross sections. 

As these first detailed calculations including relativis- 
tic and resonance effects show, the BPRM close coupling 
calculations with the large number of levels and channels 
in the intermediate energy range should considerably en- 
hance the photoionization cross section of Fe XV. Thus, 
calculations that neglect the extensive inner-shell exci- 
tation resonances are likely to underestimate the cross 
sections by large factors (as has been shown, for ex- 
ample, in the case of photoionization of neutral iron 



[10]). Although likely to be much more extensive, the 
present work emphasizes the urgent need for more elabo- 
rate BPRM calculations leading up to the L- and K-shell 
ionization thresholds, as these data would be essential re- 
quirements for photoionization modeling of astrophysical 
X-ray sources observed by CXO, XMM, and Astro-E. 
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